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ABSTRACT 

We present a measurement of the temperature-polarization angular cross power spectrum, {TE), 
of the Cosmic Microwave Background. The result is based on ^ 200 hours of data from 8 polariza- 
tion sensitive bolometers operating at 145 GHz during the 2003 flight of Boomerang. We detect a 
significant {TE) correlation in the ^-range between 50 and 950 with a statistical significance > 3.5 a. 
Contamination by polarized foreground emission and systematic effects are negligible in comparison 
with statistical uncertainty. The spectrum is consistent with previous detections and with the "con- 
cordance model" that assumes adiabatic initial conditions. This is the first measurement of {TE) 
using bolometric detectors. 

Subject headings: cosmology: cosmic microwave background 
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1. INTRODUCTION 

Since Rees pioneering work ijR.eesI 1968), polarization 
of the Cosmic Microwave Background (cmb) has been at 
the center of several theoretical studies. 

Detailed numerical predictions have been made in 
the framework of standard inflationary models with 
primordi al adiabatic and scale- i nvariant fluctuations 
(see e.g. iBond fc Efstathioulll984t iSeliak fc Zaldarriagal 
I1996D . CMB polarization data is highly useful for 
cosmology since it can shed light, for example, on 
the proce ss of reionization of the intergalacti c medium 
(see e.g. lEfstathioul Il988t iZaldarriagal Il997j) . on the 
amplitud e of the inflati o nary gravity waves back- 



ground llCrittenden et al! 119951: IZaldarriaga fc SeliakI 
Il997t IKamionkowski et al Il997|) and on the nature of 
primordial perturbations (see e.g. lSoergel & Za ldarri agal 
^97; Buchcr et al. 200lj). Moreover, polarization can 
provide further evidence for coherent acoustic oscilla- 
tions in the early universe, since in this case peaks in 
the temperature and polarization power spectra are ex- 



pected to be 180 degrees out of phase l|Kosowskvlllf 
Unfortunately, given the small amplitude of the signal, 
current cmb polarization data, while providing an im- 
portant confirmation of the standard scenario, are un- 
able to provide useful constraints on the pararn eters of 
the model. As first suggested in iCoulson et al.l (p.994), 
measuring the temperature-polarization cross correlation 
is easier, since the signal is higher and carries most of 
the cosmological information present in the polarization 
data. 

Previous detections of the temperature-polarization 
angular cross power spectrum have been obtained by 
the WMA P satelli te (Kogut et al. 2qQ3Land by the 
DASI (Ko vac et alj 2002; Leitc h et' alJ l2004l) and CBI 
l|Readhead et alJl2004,) interferometers. 

In this paper we present new observations of the 
temperature-polarization cross power spectrum of the 
Cosmic Microwave Background anisotropy obtained by 
the Boomerang experiment flown in Jan. 2003 
(BOS). Results on the temperature and polariza- 
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tion power sp ectra alone are presented in two com - 
panion papers l|,Tones et all 120051: iMontrov et~ani2005D : 

the instrument and the analysis pipeline producing the 
m aps of temperatu re and polarization are described 
in lMasi et al." f2GG 5D; the co smological parameter extrac- 
tion is in lMacTavish et all (J2()Q^). 

In the present paper we will f ollow the nota- 
tion o f 'Zald arriaga Sz Seliald l)1997j) (but see also 
IKamion kowski ^t'ldTIT997() in which polarization is ex- 
pressed as two linear combinations of spin±2 multipole 
moments which have opposite parities, the so-called E 
(electric) and B (magnetic) modes. In standard cos- 
mological models, the magnetic-type parity combination 
does not cross-correlate with temperature or the electric- 
type parity combination. The cosmological information 
in the polarization-temperature correlation is therefore 
present only in the (TE) angular power spectrum. In 
this paper we will also present constraints on (TB) as 
useful check for systematics and foregrounds. Non zero 
(TB) may also appear in exotic theories due for example 
to the presence of helical flow s in the prim ordial plasma 
at the time of recombination jPoeosian et al. 2002). 

We have performed the analysis of the BOS data using 
two completely independent pipelines, with different pro- 
cedures for the pointing solution, data cleaning, decon- 
volution, map-making and noise estimation, and differ- 
ent estimation of the calibration factors, beams, receivers 
transfer functions, polarization efficiencies and polarizer 
angles. One pipeline was developed in Italy (IT), the 
other in North America (NA). Pipelines details are in 
|Masi et al. ( 2005) and will be described furthermore in 
subsequent papers. The most important result from this 
splitting is the overall agreement, which enhances confi- 
dence in the result. A comparison of the result from the 
two pipelines allows a measure of the sensitivity of the 
result to details of the analysis. 

2. {TE) ESTIMATION 

We use data from 8 channels at 145 GHz, composed 
of 4 pairs of Polarization Sensitive Bolometers (PSB)- 
(Wl, W2), (XI, X2), (Y1,Y2) and (Zl, Z2), with effec- 
tive angular resolution of 11.5 arcminutes (full width half 
maximum, including pointing jitter). Performance an d 
characteristics of those devices are in iMasi et alJ(|2005|) . 
together with the full description of the instrument and 
of the temperature and polarization maps that are used 
for the analysis presented here. 

With polarization sensitive bolometers Boomerang 
produces maps of the three Stokes parameters, /, that 
describe fluctuations in the brightness of the radiation, 
Q and U that describe the linear polarization. The in- 
tensity of the CMB is conveniently described in terms of 
temperature fluctuations AT of a black-body respect to 
a 2.725 Kelvin black-body, and can be decomposed in 
spherical harmonics as AT{h) = J2im'^Jm^i-m{n)- Simi- 
larly the linear polarization Q + iU is decomposed using 
the spin-2 weighted basis ±2^^m 

(Q ± tU){h) - J2 («^™ T lafj ±2YUn) (1) 

irn 

thus defining the scalar field E{h) = '^gm^fm^imin) 
and the pseudo-scalar B{n) = J2em'^f-mX^ni{fi)- In the 
hypothesis that those quantities are Gaussian distributed 



and that the early Universe is isotropic, the cosmological 
information is encoded in the standard deviations and 
correlations of the coefficients: 



(XY) = {ai 



X* Y 



(2) 



where the pairs (XY) can be (TT), {EE), (BB), (TE), 
(TB) and {EB). Given the isotropy, these power spectra 
can be estimated by averaging over m at each multipole 
number £. 

Both IT and NA power spectra e stimation pip elines 
are based on the MASTER method ijHivon et alJl2002(l 
that computes the pseudo-a^m on a fraction of the sphere 
defined by the function W{n) that takes into account 
weighting and sky coverage. This yields the definition 
of mode-mode coupling kernels that depend only on the 
weighted scheme. Using an appropriate ^-binning it is 
possible to solve for the underlying angular power spec- 
tra, taking into account the binning operator, the an- 
gular resolution of the instrument, the pixelization, and 
the filtering of time stream. The quantity that is nor- 
mally used for the binning is the fiattcncd power spec- 
trum Ci = £{£ + \)Ci/2'!T. For a set of n bins indexed 

by 6, with boundaries < t^l^^^ < 4^^^\ the binning 
operator is defined as 



pb 

-III 







otherwise. 



low 



(3) 



and the power in each bin (hereafter band powers) are 
Cb = PiCg. 

The method is based on Monte Carlo simulations of 
signal-only time-streams, from simulations of the CMB, 
and of noise-only time-streams, from simulations of the 
instrument. Both simulated data-streams are processed 
in the same way as the real data, in order to take into 
account the overall effect of data filtering and partial sky 
coverage, and to estimate the noise bias to be removed 
in the power spectra estimation. 

The signal simulations are obtained from random real- 
ization of the CMB sky, in temperature and polarization, 
given an underlying cosmological model, projected in a 
time stream according to the BOOMERANG pointing so- 
lution. The noise simulations are obtained from random 
realization of the noise power spectrum, iteratively es- 
timated, taking into account noise correlation between 
channels as described in M asi et al. (2005j. 

The covariance matrix Mbb', that defines the uncer- 
tainties in the Cb determination, is estimat ed by Monte 
Carlo simulations of signal plus noise as in iHivon et alJ 
(2002). An approximation of the diagonal part of this 
matrix is given by 



'^TE.b 



(24 + l)/skyA^ 



(4) 



--TE,b 



T,b 



T,b 



Bb' 



Mi 



E.b 



E.b 



Bb' 



where /sky is the effective observed fraction of sky, Ct,!,, 
CE.b and CTE.b are the band powers of the temperature, 
the polarization and the polarization-temperature corre- 
lation respectively, TVr.b and Me^ are the band powers 
of the noise in the temperature map and in the E map, 
Bf) is the spherical harmonic transform of the beam, and 
is the bin width. 
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Fig. 1. — The {TE) power spectrum band powers for the NA 
(filled circles) and IT (open circles) pipelines. The upper part of 
the plot reports data with errorbars, the fiducial model (ACDM 
model fit to WMAP (year 1), Acbar, and CBI) as a black curve 
and the binned fiducial model as histogram. The middle and bot- 
tom plots are the results of two different consistency tests, obtained 
splitting the data in channels (WX-YZ) and in time (half 2 - half 
I) respectively. In the low-£ part of the plot is evident the ef- 
fect of a different weighting scheme between IT and NA, while at 
large multipoles the result is dominated by the same instrumental 
performances. 



2.1. Weighting 

The last flight of Boomerang, described in lMasi et all 
l|2005l). was split into three parts: a deep observation over 
0.22% of the sky, centered at RA = 82.5°, DEC = -45° 
(hereafter deep region), a shallow observation on a re- 
gion of covering the 1.8% of the sky (including the deep 
region) centered in the same coordinates (hereafter shal- 
low region), and observation of the Galactic plane that is 
not used in these power spectra analysis. Wide coverage 
and deep integration are both important for the quality 
of the result. The wide coverage of the shallow region is 
useful to reduce sample variance, the deep integration of 
the deep region to reduces the statistical noise. 

The two pipelines use different methods to combine 
the data to obtain a compromise of sample variance and 
noise. In the NA pipeline we perform independent analy- 
sis of the shallow scans and of the deep scans, computing 
the respective aj^ and af^. We then estimate four {TE) 
cross spectra, {T^Es), {TsEa), {TaEs), {TaEd}, with the 
relative correlation matrices and combine the spectra ap- 
propriately (C.R. Contaldi, in preparation). In the IT 
pipeline, a single map with all the scans is used. The 
aj^ are computed on the shallow region, the af^ on the 
deep. The effect of such a double coverage is taken into 
account in the transfer function and kernel used to derive 
the {TE) spectrum. 

2.2. Result and significance 

The results for the (TE) and {TB) power spectra are 
reported in the upper panel of Figures ^ and |2] respec- 
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Fig. 2. — The (TB) power spectrum band powers for the NA 
(filled circles) and IT (open circles) pipelines. The upper part 
of the plot reports the {TB) data with error-bars. The middle 
and bottom plots are the results of two different consistency tests, 
obtained splitting the data in channels (WX-YZ) and in time (half 
2 - half 1) respectively. 
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Fig. 3. — Likelihood of the parameters a and A£. Parameter a is 
defined as the amplitude of the (TE) fiducial model, A£ is the shift 
in multipole £ applied to the {TE) fiducial model. The continuum 
lines are for the Boomerang {TE) data, the dashed lines are for 
(TB) data, compared to the same {TE) fiducial model. In the 
central plot is reported the two dimensional likelihood, £(A£, a); 
the contours are 1, 2, and 3 tr, corresponding to 68.3, 95.4 and 
99.7% of probability. The X symbol is the expected value for (TE) 
given the fiducial model, the -f symbol is the expected value for 
(TB). In the right plot is reported the C{a) marginalizing over A£, 
and in the top plot, the C{A£) marginalizing over a. In the right 
and top plot, the gray lines are 1, 2, and 3 cr boundaries for the 
(TE) data. The {TB) data likelihood is used to test the presence 
of foregrounds and systematic effects that would affect (TE) and 
(TB) in the same way. This plot is obtained using the IT dataset 
with a binning width of 50 multipole numbers. The NA dataset 
gives a similar result. 
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TABLE 1 

BOOMERANG-03 TEMPERATURE-POLARIZATION CROSS POWER, SPECTRA BAND POWERS. 



NA IT 
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r-TE 




rTB 




rTE 




'^b 




51 


22 


32 


2 


33 


-9 


51 


2 


51 


150 


-51 


27 


-18 


27 


-71 


39 


-68 


37 


250 


40 


32 


-9 


31 


125 


46 


-73 


42 


350 


58 


28 


-16 


27 


63 


34 


29 


30 


450 


-90 


29 


7 


28 


-69 


40 


-17 


36 


550 


40 


39 


-2 


35 


20 


44 


66 


39 


650 


-18 


45 


-2 


42 


-8 


65 


91 


60 


750 


-86 


60 


-88 


58 


-88 


56 


65 


51 


850 


62 


74 


74 


72 


77 


86 


69 


82 


950 


-61 


90 


-70 


88 


-115 


72 


199 


65 


1500 


48 


81 


-133 


81 


90 


105 


-12 


100 



Units are fiK^j^j^ . The uncertainties are given by the square root of the diagonal part of the covariance matrices. The first and the last 
bins must be excluded from any analysis since can be contaminated by instrumental effects. Complete results, including window functions 
and covariance matrices, are available at http://oberon.romal.infn.it/boomcrang/b2kand http://cmb.phys.cwru.edu/boomerang. 



TABLE 2 





(TE) AND 


(TB) STATISTICS 








Test 


d.o.f. 




A 


PTE 


"^reject. 


(TE) compared to fiducial 

(TE) compared to zero 

(TB) compared to zero 

(TB) compared to fiducial (TE) 


9 
9 
9 
9 


5.81 
23.1 
4.94 
24.9 


0.48 
8.23 
0.08 
10.5 


0.62 
2.6 X 10--* 

0.92 
2.9 X 10-5 


3.5 
4.0 


(TE) compared to fiducial 
jrp (TE) compared to zero 
(TB) compared to zero 
(TB) compared to fiducial (TE) 


9 
9 
9 
9 


4.59 
20.8 
15.1 
34.5 


I. 83 

II. 3 
1.86 
20.7 


0.26 
1.2 X 10-^5 

0.16 
1.0 X 10"^ 


4.2 
6.0 


Note. — Significance of the (TE) and (TB) results respect to models. The 
contaminated by instrumental effects. 


first bin is not used 


in this analysis 


since it can be 



lively and in Tabled 

Both pipelines assume a flat shape for the power in 
each band (see equation 13)1 . For comparison with model 
band powers, the IT pipeline assumes flat band power 
window functions while the NA pipeline computes from 
the Xfaster fisher matrix estimator the band power 
window functions W^, that are used, in place of P|, to 
convert a model power spectra °^ = £{£ + 1)C^°'^/2tt 
into the theoretical band powers as 

where X [fi] = ^^'^ logarithmic integral de- 

fined in lBond et al.l l)2000() . Band powers, covariance ma- 
trices and window functions are available at the Boom- 
erang web-pages^. 

To quantify the agreement of the detection with stan- 
dard cosmology, we compare the result to C™°'^, the 
theoretical band po wers of a fi ducial model given by 
the ACDM model of iSoergel et' al. ( 2003) fit to WMAP 
(year 1), Acbar, and CBI, which we scale by a factor a 
and shift by A£. We compute the two dimensional like- 
lihood as 

C{a,Ai)(x (6) 

^ http : //oberon . romal . inf n . it/boomerang/b2k 
http : //cmb . phy s . cwru . edu/boomer ang 




where Cb',Ai are the band powers after shifting by A£ the 
power spectrum. Given the fact that the (TE) power 
spectrum crosses the zero several times, to improve the 
detection we used in this analysis a binning width of 
50 multipole numbers, and the corresponding covariance 
matrix. 

The result is reported in Figure |31 together with the 
one-dimensional likelihoods obtained by marginalization. 
For the (TE) data, the likelihood defined as above 
favours a multipole shift in the range —46 < Ai < 3 
(—38 < A£ < 20 for NA) and an amplitude in the range 
0.54 < a < 1.38 (0.40 < a < 1.30 for NA) at 95% of 
probability. The data are thus in agreement with the am- 
plitude and phase of the {TE) power spectrum predicted 
from the (TT) power spectrum under the hypothesis of 
adiabatic initial perturbations. For the (TB) data, the 
likelihood (with respect to a (TE) fiducial model) does 
not constrain the multipole shift and gives an amplitude 
in the -0.67 < a < 0.13 (-0.43 < a < 0.45 for NA) 
range at 95% of probability. 

To compare our data to a model T-Cq characterized by 
a parameter set p, we define the goodness of fit of the 
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CM 




Fig. 4. — Dust contamination. Filled circles are {TdustEBOs) ^ 
open circles (TdustBgos)- The dust contamination to {TE) is two 
order of magnitude lower than the measured {TE). The BOS data 
are from the IT pipeline. 



model as 



A ^ \ 1 ( -^(pml) 



(7) 



where £(pml) and -C(p) are the values of the likelihood 
at the maximum and for the parameters p of the model. 
In the approximation that the likelihood function £(p) 
is multivariate Gaussian near its peak, the goodness of 
fit reduces to A = and the probability of total 

exclusion is defined by the incomplete Gamma function 



PTE{K) 



1 



r(A^/2) A 



(8) 



where N is the number of parameters in the model, which 
is 2, a and in out case. A set of tests performed us- 
ing 9 bins between I = 100 and t = 1000 is reported in 
Table|21 In that Table, a represents the number of stan- 
dard deviations of a Normal distribution to have the same 
PTE as A does. The {TE) = model is rejected at 3.5 a 
(4.2 cr for IT) and the {TB) = {TE)*^ is rejected at 4.1 a 
(6.0 a for IT). The complete results of consistency with 
cosmological models and paranieter ex traction treatment 
is reported in lMacTavish et all l)2005f) . 



3. CONTROL OF FOREGROUNDS 

Polarization generated by foregrounds presents no 
global symmetry and thus is expected to contaminate 
both E and B components of the CMB in a similar 
way (see e.g. iTegm ark et al.l l2000t iTucci et all l200l 
iBacciealupj et al.l l20011. The {TB) power spectrum can 
be used to test such a contamination. The {TB)—0 re- 
sult presented in Table |21 is the main evidence that the 
{TE) result is not contaminated. 

Moreover we can directly test the contamination due 
to dust by a correlation of our data with a dust map. If 
we assume that the temperature and polarization seen 
by Boomerang are a superposition of CMB and dust, 
Tbo3 = TcMB + Tdust and £'so3 — Ecmb + Edust (and 
the same for B), then 

{TE)bo3 = {TcmbEcmb) + {TdustEdust) (9) 

where we assume that dust and CMB are not correlated 
{{TcMBEdust) = {TdustEcMB) = 0). Under the same 



TABLE 3 

Consistency tests 





Test 


d.o.f. 




P> 


NA 


{TE) temporal 
{TE) channels 
{TB) temporal 
{TB) channels 


9 
9 
9 
9 


16.1 
10.6 
13.1 
12.3 


0.065 
0.30 
0.16 
0.20 


IT 


{TE) temporal 
{TE) channels 
{TE) temporal 
{TB) channels 


9 
9 
9 
9 


14.8 
2.42 
7.55 
13.5 


0.10 
0.98 
0.58 
0.14 


Note. 


— The first bin is 


not used in 


this analysis 


since it can 



be contaminated by instrumental effects. 



assumption, the contaminating term {TdustEdust) can be 
estimated by 



{TdustEdust) — {TdustEBOz) 



(10) 



We estimate Tdust by using IRAS maps re-calibrated 
with DIRBE at 1 00 /iw. extrapolat e d to o ur wavelength 
with model 8 in Finkb einer et al.l (11999') as described 
m Masi et all (2001'). We resample the extrapolated dust 
map with the Boomerang scan strategy, and then recre- 
ate a dust map with the same time domain filtering, 
flagging, and map-making algorithm as the Boomerang 
map. As shown in Figure 0] this contaminant is compat- 
ible with zero and two orders of magnitude lower than 
the detected {TE). 

4. CONTROL OF SYSTEMATIC EFFECTS 

The standard test to detect systematic effects consists 
of splitting the data in two subsets (jackknife), making 
a differenced map using the two subsets and calculating 
the power spectra of the differenced map divided by two 
to maintain the same noise statistics as an average map. 
The result must be consistent with zero. This test is 
particularly effective because the sample variance goes 
to zero and the noise in equation Q reduces to 



(24 + l)/skyA^ 



(11) 



We performed two such tests, splitting the data in time 
and in channels. As our temporal splitting we take the 
first half of the scans on the shallow region and the first 
half of the scans on the deep region versus the second 
half of shallow plus second half of deep. As our channel 
splitting we take two PSB pairs (Wl, W2, XI, X2) versus 
the other two PSB pairs (Yl, Y2 , Zl, Z2 for focal plane 
description see iMasi et al."2005'). Results are reported 
in the bottom panels of Figures ^ and [5] and in Table 13 
showing remarkable consistency with zero. This, along 
with the {TB)=0 result presented above, gives strong 
evidence that the dataset is free from significant system- 
atics. 

4.1. Propagation of instrumental uncertainties 

Additionally, we have modelled the potential system- 
atic effects from mis-estimation of various instrumental 
characteristics using Monte Carlo simulations of signal- 
only time ordered data, processed varying those param- 
eters randomly over their range of uncertainty with a 
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10000.0 



^ 1000.0 

0~ 1 00.0 
< 



• diag(Cbb)'/2 
ir relotive colibrotion 
A transfer function 
T polarization onqle 
■ polorizotion eff. 




\ 10.0 
+ 

~ 1.0 
0.1 

200 400 600 800 1000 1200 1400 
I 

Fig. 5. — Propagation of instrumental uncertainties in the {TE) 
error-bars. The dots are the square root of the diagonal part of the 
covariance matrix, relative calibration is varied by ±0.8%, polar- 
ization efficiency by ±0.03, time constants of the transfer function 
by ±10%, the angles of the polarizers respect to the telescope frame 
by ±2°, and the beam (plotted as horizontal thicks) by ±0.3'. The 
error- bars {AC() generated by uncertainties in instrumental char- 
acteristics are one order of magnitude lower than the errors due to 
noise and sampling variance (from NA pipeline here). Those error- 
base are not treated as an increased error, but rather as a system- 
atic effect which is correlated bin-to-bin, and marginalized o ver in 
the parameter estimation as described in IBridle et al.l 



does not support models based on ac tive perturbations 
like topological defects l|Turoklll996H . In active mod- 
els the perturbations are continuously produced by the 
causal field and lead to a positive correlation in the (TE) 
spectrum. While cosmic string and textures models are 
already ruled out by the presence of peaks in the CMB 



200 




600 800 1000 1200 1400 
t 

Fig. 6. — Collection of results (TE) power spectrum from 
recent experiments. The BOOMERANG data are NA and IT, 
with alternate binnin g. WMAP dat a are from IKogu^^al 
f2003), CBI fromPIeadhead et al.' ('200f), DASI from ' LeitchTTaT 
(; 2004). PiqUE fro m de Oliveira-Co sta et al.. 1.20031D . POLAR 



from lde Oliveira-Costa et alJ 120n3aD . 



Gaussian distribution. The parameters that have been 
changed are, the relative calibration between channels 
(±0.8%), the polarization efficiency (±0.03), the bolome- 
ter time constants (±10%), the beam (±0.3'), and the an- 
gles of the polarizers axes respect to the telescope frame 
(±2°). These ranges are the uncertainties on those in- 
strumental parameters as described in Masi ct al. ( 2005) . 

As shown in Figure \S\ the potential errors from mis- 
estimation of these instrumental parameters are all at 
least one order of magnitude lower than the statistical 
error-bars of the datasct. The simulation uses the ACDM 
model of Spcrgcl ct al.. (,2003,1 fit to WMAP (year 1), 
Acbar, and CBI. 

5. CONCLUSION 

We have detected the presence of polarization of the 
CMB with high statistical significance (3.5 a combining 
the bins). This detection of (TE) con firms and improves 
previous lower fre quency detections l|Kogut et alJl2003 
iLeitch et al.ll2004l) using a completely independent tech- 
nology. The robustness of these results against fore- 
ground contamination effects is thus strengthened, and 
its cosmological origin confirmed. A summary of all mea- 
surements of the (TE) spectrum is shown in Figure |3 

The B03 (TE) data show a 2 ct anti-correlation at 
large angular scales (50 < £ < 150). This, as previ- 
ously detected by the WMAP experiment (Pcir is et al.l 
12003(1 ■ is consistent w ith the presence of superho rizon 
adiabatic fluctuations ijSoergel fc Zaldarriagalll99W) and 



temperature pow er spectrum, active models may be con- 
structed (see e.g. IDurrer 'eral]l200H) to mimic the (TT) 
data but not the {TE) spectrum. 
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